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ABSTRACT 



Values of four parameterization constants ror cne duik, 
second-order closure turbulence model of Garwooo (1976) are 
determined using ocean data from the Kortn Paciric experi- 
ment (NORPAX) Trans-Pacific ship of opportunity program 
(TRANSPAC) and atmospneric forcing from the Fleet Numerical 
Oceanography Center (FNOC) mooei. Tne annual cycle or sea 
surface temperature is simulated for eacn of tne years 
1976-1978. The difference oetween tne simulated ano tne an- 
alyzed monthly change of sea surface temperature is mini- 
mized to determine the optimal values of the mooei con- 
stants. The model is relatively insensitive to cnanges in 
the constants for this set of data. However simulations us- 
ing the values determined in tnis stuoy yieio results whicn 
are quite different from those obtained using vaiues for tne 
model constants determined from earlier studies witn differ- 
ent initialization and forcing data. 

Three statistical techniques are useo to determine tne 
sensitivity of the model to cnanges of tne model constants. 
Two techniques, bootstrapping and pooled permutation proce- 
dure, yield similar conclusions. The tniro, a standard stu- 
dent's t-test, proved to oe overly sensitive. This leaos to 
the conclusion that Gaussian statistics are not appropriate 
to studies, such as this one, for whicn tne degrees of free- 
dom are not well defined. 

The formulation rot tne adsorption or solar radiation 
with depth, in the Garwooo mooei, is a simplification of a 
two exponent model. It adequately reproduces the cnaracter- 
istics of more complex absorption models. Optimization of 
the constants of the absorption model improves tne simula- 
tions of the summer mixed layer ana the annual cycle of tne 
mixed layer. 
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1 . INTRODUCTION 

This paper aadresses three incerrelatea topics wnicn 
are fundamental to simulating the theriaooynamics of cne up- 
per ocean on aiurnai to seasonal time scales. The rirst 
topic is the determination the moael constants ror tne sec- 
ond-order closure bulk moael of turoulent oceanic Douiioary 
layers aevelopea by Garwood (1976, 1977). Tne constants are 
determined for forcing aerivea from the Fleet Nuraerical Oce- 
anographic Center (FNOC) atmospneric primitive equation moa- 
el and ocean thermal structure information oocaineu rrom tne 
Nortn Pacific Experiment (KORPAX) Trans-Pacific snip of op- 
portunity XBT (TRAMSPAC) analysis. Seconaiy, tnree statis- 
tical tests are usea to examine the variations in raooel 
skill as a function of the parameterization constants ana 
the solar radiation absorption moael. Finally, tne model 
for the absorption of solar raaiation in the upper ocean, 
which is used in the Garwood moael, is examinea ana comparea 
v/ith several other moaels. 

Turbulence closure is necessitatea oy tne process or 
Reynola's averaging wnicn generates correlations (moments) 
of tne fluctuating variaoles, resulting in more unKnowns 
than equations. The proolem can not oe resoiveu at a runoa- 
raental level with additional equations since each new equa- 
tion for a correlation introduces aaaitional nigner oraer 
moments. The system of equations must oe closea oy speciry- 
ing the highest oraer of the terras to oe retainea ana param- 
eterizing these correlations in terms of lower oraer moments 
and mean values. One may tnen characterize the various tur- 
bulence models accoraing to tne oraer at wnicn tnis closure 
is maae ana the particular parameterizations usea. Tne Gar- 
wood (1976,1977) moael discussed in the paper employs a sec- 
ond-order closure scheme. Thus, the seconu-orcier moments 
(triple correlations) are parameter izea in terras of first- 
order moments (autocorrelations ana cross correlations) ana 
zero-order moments (mean values) . 

The parameterization of the nignest order moments gen- 
erates constants of proportionality between tne nignest or- 
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der moments and functions of the lower oraer moments. Tnese 
constants must be determined experimentally before txie raooex 
can be used to simulate accurately the fluid flows or inter- 
est. The goal in turbulence modeling is to formulate tne 
parameterizations such that they are "universal” or applica- 
ble to a wide range of geophysical flows. Unfortunately, no 
parameterization can be completely universal. Parameteriza- 
tions deriveo for bounaary layer flows, for example, can not 
be expected to work well for free mixing layers, jets or 
wakes, since the physics ana tne Dounaary concitioris are 
quite different. Further, for a given parameterization, tne 
optimal values of the constants may also be expectea to oir- 
fer for different oceanic and atmospheric regimes since all 
possible pnysical processes are not incluaed in the raoael. 
In addition, the space or time resolution of the aata useo 
may not be adequate. The constants for this moael nave pre- 
viously been estimated form laboratory aata (Garwooa, 1976) 
and ocean weather ship observations at Ocean Station Papa 
14^W) (Garwood ana Aaamec, 1982) . Although these ear- 
lier values were not optimal determinations, tney will oe 
compared with the reults ootainea nere using FhOC ana 
TRANSPAC aata. 

The absorption of solar radiation in the upper ocean 
significantly influences the aiurnal cycle of the raixea 
layer, the summer minimum mixea layer deptn ana associated 
temperature maximum and the spring transition. The physics 
of the absorption of solar radiation is v/ell known; nowever, 
the distribution of yellow matter and particulate matter, 
which significantly affects tne process, is not. Further, 
an elaborate bandwidth dependent model of solar raaiation 
absorption is very computationally expensive. Thus it is 
necessary to develop an accurate out simplifiea parameteri- 
zation of the absorption process. The formulation usea in 
this model will be comparea with other methoas, incluaing an 
empirical model (Zaneveld ana Spinrad, 1980) wnich oest rep- 
resents the limited data sets that are availaoie. 
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The choice ana application or statistics wnicn are ap- 
propriate to this type of stuay is a non-criviai aspect or 
the work. Three relatively simple statistics nave oeen cno- 
sen: one Gaussian statistic ana two non-Gaussian statistics 
which are capable of generating their own prooaDiiity uis- 
tribution functions (pdf's). It is hopeo that tnis latter 
approach will reduce the problems causeo oy using sample 
populations which do not necessarily nave Gaussian pof's ana 
whose members are not indepenaent. 

The Garwooa moael is reviewed in section 2. Tne ocean 
thermal structure, surface forcing ana tne aetails of tne 
numerical simulations for this stuay are aiscussea in sec- 
tion 3. The results of the numerical experiments are aa- 
dressed in section 4 and a summary is proviaea in section 5. 



2. DESCRIPTION OF MODEL 



Since the moael has been describee in aetail oy Garwooa 
(1976 and 1977) , we will only briefly review the equations 
to indicate the moael constants which must oe aeterminea. 

Garwooa integrated the equations of motion ana the 
buoyancy equation through the mixea layer to yiela 



^<U> 

hx - f h<V> 

dt 



-uV(O) + uw'(-h) 



( 2 . 1 ) 



h 



^<V> 

at 



+ f h<U> = -vw'(O) + vw(-h) 



( 2 . 2 ) 



a<B> -7-7- -7— r ocg f 

h- < ■ = -bw(0) + bw(-h) - 



Q,- — }■ j . (2.3) 

L ^ n+6 J_ 5 j 



The vertical mean value of any variaole, X, is uerinea as 



<X> = 






X ( Z) QZ , 



where h is the mixed layer aeptn ana z is zero at the sea 
surface and positive upward. The bounaary conaitions at z = 
0 and z = -h are 



-uw'(O) = , 



oh 



uV(-h) = — 

uc 



- J - 



f 



-vw‘(0) = ^ , 

P- 

37(0) = 

f- f. 



— — . un 

vw‘(-h) = -_A_ Av — , 

at 



“TT a /» 

Dw(-h) = -_/\_ A ^ — • 

at 



Tjj and are the x ana y components of the surface stress. 

is the net surface heat flux (total heat flux, less 
solar radiation, Q^) . S is the surface salinity ana (E-P) 
is the evaporation less precipitation rate. _/\_f tne Heav- 
iside step function, is aefinea as 



_A_ = 0 for — < 0 

dt 



A 

ana _/\_ = 1 tor — > 0 . 

at 



The changes in mean velocity components ana ouoyancy tnrougn 
the entrainment zone are 

Au = U(-h) - U(-h - <T) , 

Av = V(-h) - V(-h - , 

Ab = B(-h) - B(-h - S) . 

Au . A v and Ab are the cnanges in the mean properties 

across the entrainment zone ana ah/at is the race of raixea 
layer deepening. 

The turbulence equations integrated chrougn the mixea 
layer are 



d(h<u^ + v^> 
dt 2 



= ra_u5 + 

3 * 



■(AU)^ + (AV)'^'lah 






— • - m.,<e>'A(<e> - 3<w^>) 



at 



2m, r — , m< 

’ <e>^ + -^fh 

m. 



(2.4) 



<e> 



d(h<w2»>) h bV(0) 



at 



- abS + Sii-ff'iki . f 

at a.C,U,L 2 4 



Q- ( z ' ) az ' I aj 



m. 



+ m.,<e>2i(<e> - 3<w3»>) - — 

* 3 L 



— 2 % 1 — 

<e>^ + -^f n <e> . 



J i 
(2.5) 



m. 



Integrating the total turbulent kinetic energy (TKE) tnrougn 
the entrainment zone provides the final equation neeaea to 
close the system ^of equations, 
dh 2ra^(w lri<T> 

dt hAB + 2<^> 

<e> is the vertically averagea TKE 



(2.6) 
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<e> 



<uJ»> + <w^> + <w2*> 



2 
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The left hand siae of (2.4) is the local rate of cnanoe of 
the vertically averaged horizontal TKE in tne mixea layer. 

The time scale for the terms on the left nano siae of 
(2.3) and (2.4) is approximately 100 seconds for tne oceanic 
planetary boundary layer under normal circumstances (de 
Szoeke and Rhines, 1976). Since we are interestea in montn- 
ly ana longer time scales in tnis stuay, tnese terras will oe 
neglected. The first term on the rignt nano siae of (2.3) 
is a parameterization of the surface wina stress acting on 
the shear of the mean current, the proauction aue to oreax- 
ing surface waves and the transport of norizontal TKE. u^^ 
is the ocean friction velocity, given Dy = jTl/p^ana 
used as a scale for the rms velocity fluctuations. Tne sec- 
ond term on the right hand side of (2.3) is tne momentum 
flux associatea with the entrainment of nonturouient fluia 
at the base of the mixed layer. The momentum flux aue to 
internal waves is neglectea. The thira term is a parameter- 
ization of the pressure rate-of-strain terms wnicn act to 
reduce the anisotropy of the turbulent flow (Rotta,1951; 
Lumley and Khajeh-Nouri,1974) . Tne last term is a parame- 
terization of the aissipation rate. The aissipation rate 
can be estimatea from the rate at whicn tne large scale ea- 
dies supply energy to the smaller scales. As aiscussed in 
Tennekes ana Lumley (1972) , the dissipation rate is 



where T^ is the dissipation time scale. There are two time 
scales for the oceanic mixed layer; the time scale imposea 
by planetary rotation, T^^^l/f ana the time scale rormea 
from the turbulent velocity ana the length scale of tne 
large scale turbulent flow, h/u^. The first scale is 
the characteristic time scale of vortex stretcning aue to 
the planetary rotation. The second is the time required ror 
a typical large scale eaay to overturn. A feature of Gar- 
wood's model is that both of tnese time scales are incorpo- 
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rated into the parameterization of the aissipacion time 
scale by specifying that is scaled with the smaller of 
or T^. To incorporate both time scales the aissipation time 

scale is defined as 

•f *1 “I 

Te = Tr + Tt . 

In (2.5) , the first two terms on the rignt nano siae 
are the surface and entrainment buoyancy fluxes. Tne txiird 
term is the absorption of solar radiation witn aeptn. Tne 
remaining terras are identical to those in (2.4). (2.1) 

through (2.6) form a complete set of equations for tne ae- 
pendent variables <U>, <V>, <B>, <u^ + v*«>, <w^> ana n. 

The depth-depenaent absorption of soiar rauiation is 
particularly important in shallow, summer mixea layers. In 
these cases the amount of solar raaiation absoroea in tne 
mixea layer significantly affects the mixea layer teiapera- 
ture, while the absorption below the mixea layer arrects the 
entrainment rate and, therefore, the mixea layer aeptn. Tne 
flux of solar radiation in the ocean ean oe written as 
QJz) = Q^<o)| dA 

where A is the wavelength ana W(A) is tne optical water type 
(Jerlov, 1976). Various parameter izations (see Taole 1) 
have been proposed for the aosorption of solar raaiation. 



TABLE 1 

Parameter izations of solar radiation in tne ocean 



Parameter izations 



Source 



Qj(z) = 

Qj(z) = (^c.)(r, + rj e'^»*) 

Q^z) = r; e'*i^ 

Qj(z) = - a tan”* (ox)) 



Denman, 1973 

Paulson ana Simpson, 1977 
Konao, 1979 

Zanevela and Spinraa, 1980 



Q^z) = lz| <1 Garwood, 1976 

= Q‘»)re^^ lz| >1 

Q«u)is the solar radiation at the sea surface, the 0 ' s are the 
extinction coefficients ana the r's are constants (0 < r < 1) . 



The first three of these parameterizations assume tnat a 
significant fraction of the solar radiation aosorption can 
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be modeled by using one, two or many wavelength Danos, re- 
spectively. The fourth parameterization was ootainea empir- 
ically using an extensive data set of absorptance ana trans- 
mittance measurements (Spinrad et. al. , 1979) . Tne con- 

stants in the fourth formula in Taole 1 aepena on tne water 
type. For this study the appropriate water type is Type II 
(Zaneveld, personal communication), and the constants are 
= 0.067 m/s, a = 0.4158, ana b = 3.9865 m/s. Simpson anc 
Dickey (1981) comparea the first four parameter izations in 
Table 1 ana founa that using a single exponential mouei is 
acceptable only at very high wina speeas (U,^ > 20 m/s) . 
They also founa that the empirical formulation of Zaneveia 
and Spinrad (1980) and the double exponential mooei or Simp- 
son ana Paulson (1977) produce nearly identical temperature 
profiles, and that the complexity of tne spectral aecomposi- 
tion form is not necessary in mixea layer raoaeling. 

The fifth parameterization, usea in the Garwooa moael, 
is a douDle exponential moael for whicn one of the extinc- 
tion depths approaches zero. A fixed fraction (1-r) of Q 
is absorbed in the first meter ana the remaining flux is ao- 
sorbed exponentially for |z| > Im. Values of r ana o must 
be specified. This parameterization will oe comparea witn 
Zaneveld' s empirical formulation. 

In (2.4) to (2.6), m, througn m^ are tne constants 
which must be determined from laooratory ana geophysical 
data. Of these only four are inaepenaent, ana only three (ir^ 
and p^ = m^/m, ) 

are critical to moaeling the annual cycle of the upper 
ocean mixed layer. In this stuay, we use heat ana momentum 
fluxes derived from the FNOC atmospneric fielas ana tne 
TRANSPAC ocean data over a range of latitudes ana longituaes 
to estimate the optimal values for m^, p^= m^/m^ , r ana 

3. DESCRIPTION OF NUMERICAL EXPERIMENTS 

Choosing an annual time perioa for this stuay requires 
that the moael be able to simulate four signir leant pnenome- 
na each of which requires aifferent responses from tne moa- 
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el. The fall deepening is particularly sensitive to ruj (eq. 
2.4), the parameterizacion of surface wind mixing. Tne win- 
ter minimum mixed layer temperature is sensitive to p^ («iqs. 
2. 4, 2. 5), the dissipation parameterization. The spring 
transition aepends on the interaction of m ana tne surrace 
heat fluxes. Finally, the summer maximum mixea layer temp- 
erature depenas on r ana ^ in the parameterization of solar 
raaiation absorption. 

A grid covering a large area of the Kortn Paciric rrom 
34*N to 46^N in 4^ increments ana from 170*W to 140^17 in 10^ 
increments was chosen for this stuay. This area is large 
enough to encompass various oceanic tnermai structure ana 
atmospheric forcing regimes but aoes not inciuae areas wnicn 
are strongly influenced oy eastern or western boundary cur- 
rents where horizontal processes (e.g. aavection ana aiver- 
gence) are likely to be important. The north-south extent 
of the grid incorporates only the area whicn nas consistent 
data coverage throughout the entire year. Regions to tne 
north (south) of these latitudes have fewer aata points cur- 
ing the winter (summer) aue to the latitudinal snift of ship 
tracks with the seasons. 

The specific quantity chosen to compare the moaei re- 
sults with the data is the montnly change in sea surface 
temperature (SST) averaged over an annual time perioa ana 
over the grid. The monthly change in SST is more difficult 
to simulate than the SST since the large scale nortn-soutn 
SST gradient in the nortn Pacific contains more of tne vari- 
ance than that associatea with the month-to-moxith cnange in 
the SST field. J..v. v^.-<_nti cy to ot= luinimizea can ue written 
ss 

/i w 

^ ' r I r I 

where M is the number of gria points, N is the numoer of 
months used in the average, ana A.T^ ana At^^ are the montnly 
cnanges in the simulated and analyzed SST respectively. 
This quantity was choosen ratner tnan an rms airference 
since it allows the study of some additional statisiticai 
techniques which will oe reportea elsewhere. 
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The ocean thermal structure tields whicn are useu in 
the initialization of the model ana in the computacion of 
the model error were obtaineo from the NORPAX TRANSPAC anal- 
ysis (V7hite ana Bernstein, 1979) . Tne optimal interpolation 
analysis uses the temperature ooservations during a montn to 
estimate the temperatures at each grid point. Horizontal 
analyses are prepared, independently, at 0, 20, 40, 60, 80, 
100, 120, 150, 200, 250, 300, and 400 meters. These rieids 
are interpolated to 5m depth intervals using a combination 
of linear ana exponential interpolation (Elsoerry, et. ai. , 
1982) . The model is initialized on January 15th ana inte- 
grated through December 15th for 1976, 1977 ana 1978. The 
model generated monthly SST used in (3.1) is a five aay av- 
erage centered on the 15th of each month. 

The surface boundary conditions requirea for tne Gar- 
wood mooel are the total heat flux (Cy) , the solar radiation 
(Q^) and the wind speed (W) . These variables are deriveu 
from fields generated by the FNOC atmospheric preuiction ana 
analysis system. Instantaneous values of ana are ar- 
chived at 12-h intervals and W at 6-h intervals on a 63 x 63 
polar stereographic grid of the northern hemisphere. Biii- 
near interpolation is used to map these fielas onto a 4 la- 
titude by 10* longitude grid chosen for this study. Tne 
values are interpolated to 1-n intervals (Gallacner, 1979) 
and Q.J. is corrected to remove a systematic long term bias 
(Elsberry et. al., 1982). 

The procedure used to fina the oest fit of tne moael 
simulations to the analyzea fields was to first set r = 0.5 
and if = 0.1 m"* and minimize (3.1) in (m^, p^) space ror 
each year. Next, r was nela constant ana m^ , ana p^ were 
varied to achieve the best fit to tne 1977 annual cycie. 
This value of was compared to the optimal ^ found indepen- 
dently by comparison with Zanevela's empirical formula. In 
the comparison with Zaneveid's formula botn r ana ^ are var- 
ied. 
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4. DISCUSSION OF RESULTS 

As a preliminary test, the 1977 annual cycle v/as usea 
to determine the distribution of the model error in 
space. A broad but physically reasonaole range or anu p^ 
was tested using widely spaced values of ra_, p . Tnese sim- 
ulations indicate a local minimum of E in che neignoornooo 
of m^ = 2.0, Pj = 4.0. The minima for all three years were 
located by testing the region centered on p = 4.0, m = 2.0 
(Figs, la, b and c) . For each year che minimum lies on tne 
line p^= 4.0, but the value of m^ varies oecween 1.5 ror 

1978 ana 4.6 for 1976 . 

4.1 STATISTICS QZ VARIATIONS M lii£ MODEL CONSTANTS 

The student's t-test is used to determine wnecner cne 
variations in E are stacistically significant. Tne null ny- 
pothesis, , is E(l) - E(2) = 0 wnere 1 ano 2 represenc 

any two values of the spatial and temporal average moael - 
analysis differences defined oy (3.1). The nypochesis is 
accepted if 

-t < t < t 

o(/2,m+n-2 ot^2,m+n-2 

where t is the student variate, oC is tne significance level 
and m and n are the numoer of points in the samples (m = 16, 
n =10 ; MxN = 160, for these scuoies) . Tne student variate 
is defined as 

E(l) - E(2) 
sp (1/m + 1/n) ' 

where sp, the pooled variance, is 

(m-1) si + (n-1) s2 

sp = — . 

(m+n-2) 

si and s2 are the variances of E(l) ano E(2) respectively. 

For the values in 1976 (Fig. la), tne student's t-cest 
indicates that all non-adjacent points are signir icantly 
different at the 90% confidence level (which will be used 
for all cases) . Furthermore, the only ad 3 acent pairs of (m_j 
/ Pj) points which do not differ significantly are (3.0, 4.0 
and 3.5, 4.0), (3.2, 5.0 and 3.5, 4.0) ana (3.2, 5.0 and 



10 



3.2, 4.0). In the 1977 siraulacions, tne results are .aucn 

less sensitive to variations of the parameters. Tne porno m_ 
= 1.5 was significantly cifferent trom the poinos hi^ =2.2 
ana = 2.4 v/ith p^ = 4.0. None of ohe other points were 

significantly Different for 1977. In 1978, wicn p^ = 4.0, 
the point m^ = 1.8 differs signit icanoly from ruj = 1.5, 1.6, 
ana 3.0 and m_j = 2.6 is significantly aiffereno from = 
1.4 and m^ = 1.5. If the constants are see ao p^ = 4.0, 

= 2.2, the raoael-analysis Differences are signiricant, ac 
the 90% level, among the three years. Furtnermore, one min- 
imum in the parameter space for eacn year is signif icanoly 
different from the minimum in eacn of one ooner owo years. 

The student's t-oest assumes ohao the samples are arawn 
randomly from a normal population. Thus the aata poinos are 
assumed to be indepenaent of eacn ooher. This assumpoion 
can be tested by pooling the /\t ' s from two seos of moaei 
simulations. A probability aistribution function (pat) can 
be generated by ranaomly arawing samples from tne new popu- 
lation. This pdf can then be usea to determine tne signifi- 
cance levels. Using this method aoes not require thao one 
samples have a Gaussian pdf or thao one points oe inaepen- 
dent. An example of this procedure, tne pool ana permuoa- 
tion procedure (PPP) is descrioea oy Preisenaorf er ana hoo- 
ley (1982). First, the two data secs are pooiea inoo one 
set where the first m elements are the first sample of At's 
and the following n elements are the secona sample. Tne 
next step is to permute ohe m + n aata poinos to consoruco a 
new pooled data seo. The first m elemenos oecome a new 
first sample ana the following n elements become one secono 
sample. The test variable t may ohen oe calculaoea ror 
these two new aata sets. These steps are now repeaoea for a 
total of 100 random experiments yielding 100 values of t. 
The t are then orderea in increasing value v/ion o(l) oeing 
the smallest ana t(lOO) tne largest. Then the values of 
t{5) ana t{95) now yield the 90% significance levels airecc- 
ly, which eliminates the assumptions required to use ohe 
student's t-test. 
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A similar method callea 'bootstrapping' can oe usea co 
generate the pdf ( Preisendorf er , personal communicacion) . 
In bootstrapping, the selections are maoe randomly rrom the 
entire data set, whereas only the permutations of cne pooiea 
population are usea in the PPP. Therefore, since repeti- 
tions of the data points may enter into the pooled set using 
bootstrapping, slightly broader reference aistr ioutions may 
be generated. Thus, the null nypothesis would tend to oe 
accepted more often using bootstrapping than using tne PPP. 

The critical points for accepting the null hypothesis 
may be quite different from the student's t-test if the 
original data sets do not have approximately Gaussian pafs. 
If this is the case, fewer runs will be statistically dif- 
ferent using the oootstrapping method and PPP. Results are 
taoulated in Taoles 2, 3 ana 4 for cases in whicn the stu- 
dent's t-test rejected the null nypothesis at tne 90% conri- 
dence level. Variations in ra^ ana p^ within a given year 
(Table 2) do not lead to any significant differences except 
using the PPP between m^ = 2.0 and m^ = 3.5 for p^ = 4.0 in 
1976. 



TABLE 2 



Comparison of PPP ana Bootstrapping 



Year 



1976 

1976 

1976 

1977 

1978 



(nijfPj ) 


( nij f P3 ) 


t 


(3.5,4.0; 


(2. 0,4.0) 


-11.49 


3. 2, 4.0 


3. 2, 5.0 


- 3.09 


3. 2, 4.0 


2. 2, 4.0 


- 8.11 


1.5,4 .0 


2. 4, 4.0 


1.86 


(1.5, 4.0 


(3. 0,4.0 


- 5.78 



PPP 

c(5) t(95) 



-11.05 S.66* 

•11.96 10.01 
-12.31 8.02 

- 6.91 8.97 

- 9.46 7.73 



Bootstrapping 
t(5) tTss) 



-12. a8 9.72 
-16.27 10.69 
-12.15 9.62 

- 8.58 7.93 

- 7.65 7.73 



Student variate ana 90% confidence limits as determined oy cne 
bootstrapping method ana by PPP for the cases for waicn tne scuae 
t-test rejected the null nypothesis at the 90% levei wicnin each 
particular year. The entry marked with an aster isk was tne only- 
case for which either test determined that the points were 
significantly different. 



Holding the parameters constant (p^ = 4.0, m 2 = 2.2) and 
testing between years (Table 3), the PPP test leads to the 
conclusion that the results for 1976 ana 1977 are signifi- 
cantly different, but tne oootstrapping tecnnique aid not. 
Both procedures found 1976 ana 1978 to be significantly aif- 
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4.0 out not for m 



ferent for = 3.0, = 



= 2 . 2 , 



= 4.0 



3 - t'j 

(Table 3). Results for 1977 ana 1978 are not statiscicaliy 
different. Testing for differences between minima or eacn 
year (Taole 4) , both tecnniques founo than the 1976 resuics 
are significantly different from 1977 ana 1978, out 1977 ana 
1978 are not significantly different. 



TABLE 3 

Comparison of PPP and Bootstrapping Between Years 



Years 



1976 

1976 

1977 
1976 



vs 

vs 

vs 

vs 



1977 

1978 
1978 
1978 



(m, 



f Pj ) 



2 . 2,4 . o ; 
2 . 2 , 4 .0 
2 . 2 , 4 .0 
3 . 0 , 4.0 



PPP 

t t(5) t(95) 

- 8.01 - 7.13 9.38* 

- 4.88 - 7.84 7.86 

• 1.83 -10.51 7.58 

■13.16 -12.08 10.34* 



Bootscrapoing 
t(5) t(§5) 

• 8.82 9.21 

• 7.75 7.43 

• 7.61 8.81 

■ 9.53 8.93* 



Student variate and 90% confiaence limits as aecerniined by cne 
bootstrapping methoa ana by PPP for the cases for wnich cne stuaenc's 
t-test rejected the null nypochesis at tne 90% level. Tne entries 
marked with asterisks were the only cases for wnicn either test 
determined that the points were significantly aifferent. 



TABLE 4 

Comparison of PPP ana Bootstrapping Between Minima 



PPP Bootstrapping 

Years t t(5) t(95) t(5) cT95) 



1976 


VS 


1977 


-17.49 


- 6.59 


7.47* 


-12.18 


5.86* 


1976 


vs 


1978 


- 9.55 


- 7.42 


7.44* 


- 9.22 


6.90* 


1977 


vs 


1978 


2.42 


- 7.88 


6.00 


- 9.35 


6.47 



Student variate ana 90% confidence limits as aeterminea by the 
bootstraoping method and by PPP for the minima oetween eacn or 
the three years. These (m,,p ) minima are: 1976 (4. 6,4 .0), • 

1977 (2. 4 , 4.0); ana 197M i?5 , 4 . 0 ) . 

The entries marked with asterisks were cases for wnicn the test 
determined that the points were significantly aifrerent. 



For most cases, the bootstrapping methoa yielas smaller 
values of t(5%) ana larger values of c(95%) than aoes the 
PPP. In adaition, both techniques yield smaller values of 
t(5%) and larger values of t(95%) than the scuaent's t-cest. 
Again, this implies that the bootstrapping metnoa ana tne 
PPP will fina two cases to be significantly aifferent less 
often than will the stuaent's t-test. This large aifterence 
in the critical values usea for the different types of tests 
indicates that the assumption of a normal population is not 
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valid and the data points are not indepenaenc. Tnus cne 
student's t-test is rejeccing when it shoulo De acceptea. 

4.2 EXAMINATION Q£ SOLAR RADIATION MODEL 

Two approaches were used to determine optimal values 
for the constants, r and in Garwooo's solar raaiacion pa- 
rameterization. In the first method, regression analysis 
was used co determine the values of r ana ^ cnac proviaeo 
the best agreement between Garwood's ana Zaneveld's parame- 
ter izations. Then r ana ^ were cnangea ana m^ ana ^ were 
varied to minimize E for 1977. 

Since Garwooo's parameterization is an approximacion to 
a two bandwidth model of solar raaiation aosorption, a non- 
linear regression was performed to compare Zanevela's param- 
eterization with that of Paulson ana Simpson (see Taole 1) . 
Zaneveld's model of solar raaiation absorption is an empiri- 
cal representation of the extensive aosorptance ana trans- 
mittance data set collectea by Spinrad et. al. (1979); tnus 
it is assumed that agreement with Zanevela's moael is tanca- 
raount to agreement with the data set. The values wnicn gave 
the best agreement in the upper 10m were (r =0.4, = 

0.08m ) and (r = 0.6 and ^ = 3.0m ). The assumption mat 
one of the two extinction depths is much smaller than cne 
other is verified in that is approximately cwo orders of 

magnitude larger than . The upper 10m was cnosen since 
the absorption of solar raaiation above the extinction 
depth, is most critical to the mixea layer hear ana ouoyancy 
budgets. Linear regression analysis was also usea to com- 
pare the logarithm of Garwood's parameter izacion to cne log- 
arithm of Zanevela's empirical rormulation ror -10m < z < 
-Im. That is, the equations 

In ( r ) ■¥ K z 

and 

^ z + ln(l - a tan”*(bx)) 

were compared to obtain the values of r ana ^ v/nicn gave che 
best agreement. The best agreemenc was obcainea for (r = 
0.4, ^ = 0.08m'*), wnich again verifies thac Garwooo's param- 
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eterization approximates che aouble exponential mooeir wnicn 
was found to aoequately aescrioe tne aosorption of solar ra- 
diation witn respect to raixea layer cynamics {Simpson ana 
Dickey, 1981) . 

We are concerneo with the effects of tne parameteriza- 
tion of solar radiation on the neat Duaget of tiie mixeu lay- 
er. With this in mind, Zanevela's ana Garwooa's parameteri- 
zations were compared by calculating the cnange in 
temperature of an initially homogeneous water column aue to 
the aosorption of solar raaiation. The initial temperature 
of the water column was 12*C ana the equation 

3 T 3 Qjz) 

3 z ^Cp 3z 

was integrated for 12 hours. The surface solar raaiation 
had a cosine dependence in time v/itn a peak at local noon. 
The temperature profiles after 12 hours for three cnoices of 
r and o are shown in Fig. 2a. The results for (r = 0.4, ^ = 
0.08m ) and (r = 0.45, 0 = 0.095m ) both overestimate tne 

temperature at Im, underestimate the temperature at 2m ana 
slightly overestimate tne temperature from 4 to 10m. Tne 
set (r = 0.5, = 0.11m*) estimates tne temperature at Im 

and 2m well but overestimates the temperature from 2 to 10m. 
The change in heat content is best estimatea oy ( r = 0.45, 3” 
= 0.095m*), since the uncer estimate of the temperature in 
the first 2m is approximately balancea oy tne overestimate 
of the 2 to 10m temperatures. Next, the value of r, tne 
fraction of solar raaiation wnicn penetrates 1 m, was neia 
constant ana o was variea. For r = 0.5 ana 0.45, ^ = 0.11m 
and o = 0.095m"*, respectively, gave tne uest agreement oe- 
tween Zanevela's ana Garwooa's parameterizations. In Fig. 
2b, profiles are plotted for r = 0.50 ana JT = O.lOni ', 
0.15m’*, and 0.20m'*. Fig. 2c is similar to Fig. 2D, out for 
r = 0.45. In both of these figures, all three choices of ^ 
give good agreement near the surface v/nere the graaient witn 
aepth is very large. Compared witn Zanevela's curve, tne 
curve for which ^ = 0.20ra* turns too soon, wnicn implies 
that the scale depth for absorption is too aeep. Tne curve 
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for which = O.lOnT* has approximately equal area aoove ana 
below Zaneveld's curve; which makes the integracea hearing 
in the upper 10m the same. 

The annual cycle for 1977 was simulatec with r = 0.50 
ana ^ = 0.10m"*, 0.15m * ana 0.20m * for a range of values of ra^ 
and p_ . The values of E were higher for it = 0.15m * ana V = 
0.20m"' than for if = 0.10m and the minima were aisplacea to 
higher values for higher ^ values. The cnange in m^ can 
be explained since larger ^ values imply shallower extinc- 
tion depths. Therefore, more solar raaiacion is aosorbea 
near the surface ana larger surface stresses are requireu to 
mix the more stable water column. Using tne Dootstrapping 
method or the PPP, there is no significant difference among 
results for V = 0.10m *, = 0.15m* ana = 0.20m”* ac the 90% 

confidence level, where (m^= 2.5, p^= 4.0). 

4.3 CO M P A RI SON WI.TR PREVIOUSLY C ALC .UL M £B VA L.U.BS 

The optimal values of the model conscancs found in chis 
study are quite different from the values previously used 
and emphasize the need for careful evaluation of tne con- 
stants in different conditions. The values of these con- 
stants determined by Garwooa (1976) were m_ = 7.6, p_ = m^/m, 

V -I J J 

= 1, r =0.5 and 0 = 0.1m . The Garwooa model nas previ- 
ously been tested with OS Papa data by Garwooa ana Aaamec 
(1982) , who specified m- = 0.75, p = 1.0, r = 0.5, ana if = 
0.2m . The ocean thermal structure and surface flux infor- 
mation used in the present sruay is quire uifferent from 
rhat used in the OS Papa simulations. The ocean thermal 
structure used in the OS Papa cases was aerivea from XBT oo- 
servations taken at the same site approximately every two 
days. The thermal structure for this study was oorainea 
from monthly objective analyses of XBT oDservacions on hori- 
zontal surfaces of constant depth. Also, the XBT oDserva- 
tions can resolve temperature graaients wirh ^z < 5m, 

whereas the TRANSPAC analyses ao not resolve graaients wicn 
^z < 20m. In the OS Papa simulations of Garwooa ana Aaamec 
(1982) , the surface stress and rhe hear fluxes were calcu- 
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lated from 3-n observations of SST, air cemperacure , 
point temperature, cloud cover, wind speed ana aireccion. 
The v/ina speed ana direction were aajustea to 10 m. For 
this study, the surface stress is calculated from a mouei 
producea surface wino speea adjuscea to 19 m. Tne cocal 
heat flux ano solar raciation flux are calculacea froui i.ioaei 
heat and moisture fields. 

Model sensitivicy co variations of cne values of cne 
model constants may be juuged by comparing time series of 
SST generated using Different values of the constants with a 
time series of the analyzea SSTs. Tne improvement in moael 
performance using the values aeterruinea in tnis stuuy cora- 
pared with either the Garwooa values of tne constants or tne 
OS Papa values is illustratec in Figs. 3a, d ana c. In Fig. 
3a, using values which are typical for tnis study (m^ =2.0, 
p^ = 2.0, r = 0.5, ^ = 0.1 m')f tne moael accurately simu- 

lates the annual cycle of sea surface temperature. Using 
the OS Papa values {rn_= 0.75, 1.0, r = 0.5, ^ = 

-I J 3 

0.2m ), the maximum summer temperature is too hign aue to 
too little wina stirring (m^ = 0.75) and aue to tne snaliow- 
er extinction depth = 0.2m’*) used in aosorption of tne 

solar radiation {Fig. 3b). Further mere is less variabili- 
ty in the summertime mixea layer depth ana tne rate of tne 
fall mixed layer deepening is less than in Fig. 3a. In Fig. 
3c, using the values aetermined by Garwooa (m^ =7.6, Pj “ 
1.0, r = 0.5, o= 0.1m ) , the moael aoes not warm enougn in 
the summer. This is due to tne fact that m^ is too large, 
causing the effect of the wina stirring to oe too great axia 
Keeping the mixed layer mixed too aeep. Tne variability in 
tne summertime mixed layer aepth ana tne rate of tne rail 
deepening is much greater than in Fig. 3a. Also tne winter 
mixed layer is too deep comparea to Fig. 3a. Anotner impor- 
tant difference among the three simulations is aeptn of tne 
summer mixed layer ana the time of the occurence of tne 
spring transition. The spring transition, wnicn is tne 
abrupt and severe retreat of the mixea layer, occurs near 
the end of February for tne simulations using tne constaxjts 
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determined in this study and the OS Papa constants. Hov/- 
ever , in the simulation using the Garwooa constants, the 
transition is delayed until mid-April and the upper ocean 
never achieves the shallow mixed layer ana strongly strati- 
fied seasonal thermocline that are characteristic of tne 
summer mixed layer in raid-latitudes. 



5. SUMMARY AND CONCLUSIONS 

The optimal values of the constants for any turoulence 
model are expected to vary for different specifications of 
the atmospheric forcing ana the ocean thermal structure. We 
have shown that the values of the moael constants estimated 
from laboratory data, and Ocean Weather Snip observations 
yield quite different results, when applied to the same sim- 
ulation (Figs. 3a, 3b and 3c) , than the "optimal" values de- 
termined here. Whether it would be oetter to alter tne val- 
ues of the exchange coefficients for the surface tiuxes 
rather than the values of the closure constants is largely a 
philosophical question. Knowledge of surface fluxes ana 
geophysical turbulence is not yet adequate to distinguish 
between the two approaches. We nave cnosen the second ap- 
proach largely for its inherent simplicity ano to demon- 
strate the importance of understanding the potential inter- 
actions Detween the model and the data usea to initialize 
and force the model. 

The optimal values may also cnange for different atmos- 
pheric and oceanic regimes, because all possiDle pnysicai 
processes are not parameterized in the moael. We nave snown 
that the Garwood model is relatively insensitive to varia- 
tions in the model constants within a fairiy oroad region of 
parameter space. However, the year-to-year variations in 
the optimal values of the model constants are significant. 
In particular, 1976 was found to be significantly dirrerent 
than 1977 and 1978 by all three statistical tests employed 
in this study. This reflects the fact that tne atmospneric 
forcing and the response of the ocean thermal structure were 
quite different among the three years (White, et. ai. , 1980; 
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Haney, 1980; ana Elsberry, 1983). An improvement in raoaei 
skill could be achieved by calculating the optimal values or 
the model constants for each year. Because one acmost/neric 
forcing and the ocean thermal structure are calculated in 
the same manner for all three years, the variations in tne 
optimal values of the constants probaoly represent tne er- 
fects of physics not incorporateo in the mooel anu noise in 
the analyses of ocean thermal structure ana tne surface 
fluxes. For these reasons, the same values will oe useo for 
all experiments using this type of forcing in the North Pa- 
cific. The values, p = 4.0, m_ = 2.0, r = 0.5, ana = 0.1 

-I u ^ 

m , were chosen as oeing representative or tne swarms or 
points from the three years tnat oo not cirrer signir icantly 
at the 90% confidence level. 

Three statistical tools, tne pool ana permutation pro- 
cedure, the bootstrapping technique and the stuaent's t- 
test, were used to determine wnether cifferent parameter 
values gave statistically different model results. In gen- 
eral, the bootstrapping technique generated siigntiy oroauer 
reference distributions; hov;ever, it influenced tne infer- 
ences drawn from the tests in only one case (1976; = 3.5 
vs. m_j = 2.0). The student's c-test, on tne otner nano, 
lead to rejection of the null hypotnesis in far too raany 
cases. Tnis points out tne hazard of assuming Gaussian 
pdf's for this type of data. V7e conclude tnat statistical 
tests wnich put less weight on the independence of the data 
points are more useful for studies of this nature. 

The aosorption of solar radiation witn aeptn is not 
considered one of the fundamental aspects of geopn^^sical 
turbulence modeling. However, we have shown that its proper 
specification is important for simulations of tne annual cy- 
cle and crucial to simulations of the summer mixed layer. 
The parameterization employed in the Garwood model reproduc- 
es the results of more complicated formulations ana, witn 
optimized coefficients, improves the summer anc annual simu- 
lation of the upper ocean temperature. 
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LIST OF FIGURES 



Fig. la Mean model-analysis temperature airrerence ( 
C) for 1976 computea using (3.1). Tne minimum is 
markea by an asterisK. 

Fig. lb Similar to Fig. la except for 1977. 

Fig. Ic Similar to Fig. la except for 1978. 

Fig. 2a Temperature profiles in tne upper 10m of tne 
ocean for Zanevelo's empirical formulation anc Gar- 
wood's formulation using three r.a pairs: (0.40, 
0.080), (0.45, 0.095) ana (0.50t O.lIO). For eacn 
value of r a regression was pertormea against Zane- 
vela's empirical formulation to fina tne optin.al val- 
ue of o. 

Fig. 2b Temperature profiles in the upper 10m or tne 
ocean for Zaneveld's empirical rormuiatioii ana ror 
three (r,o) pairs or Garwooo's parameterization. In 
these, r was set to 0.50 ana « v/as variea (0.10, 
0.15, ana 0.20) . 

Fig. 2c Similar to Fig. 2b except for r = .45. 

Fig. 3a Time series of aaily mixea layer temperature 
ana aepth as predicted by the model for 1977 at (3b 
M, 160*W) . The X's aenote the analyzea SST ror cnat 
month. The moael constants aeterminea in tnis stuoy 
(m^ = 2.0, p, = 4.0, r = 0.5, i = 0.1m"') were usea in 
the simulation. 

Fig. 3b Similar to Fic. 4a except using OS Papa ragoel 
constants (m, = 0.75, p^ = 1.0, r = O.o, ^ = 0.2m'). 

Fig. 3c Similar to Fig. 4a except using Garwooo 
11976) model constants (m_ =7.6, p_ =1.0, r = 0.5, 
if = 1 .Onf) . * ^ 
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Fig. 2a Temperature proriies in tne upper 10m or tne 
ocean for Zaneveio's empirical f ormu^iation ana Gar- 
wood's formulation using tnree r.e pairs: (0.40, 
0.080), (0.45, 0.095) ana (0.50^ O.llO). For eacn 
value of r a regression was oerrormea against Zane- 
veia ' s yOmpir icai formulation to rinu tne optimal val- 
ue of o. 
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2b Temperature profiles in the upper 10m of tne 
for^ Zaneveia's empirical formulation anci tor 
(r,0) pairs of Garwooa's parameterization. In 
r was set co 0.50 ana if was variea (0.10, 
ana 0.20) . 
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